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XV. 


ON CAUSTICS* 
PART FIRST 
[1824] 


To the Reverend Archdeacon Brinkley, Professor of Astronomy in the University of Dublin 
and President of the Royal Irish Academy, this first paper on Caustics is respectfully inscribed. 


The problems of Optics, considered mathematically, relate for the most part to the inter- 
sections of the rays of light, proceeding from known surfaces, according to known laws. 

In the present paper it is proposed to investigate some general properties common to all such 
Systems of Rays, and independent of the particular surface, or particular law. It is intended, 
in another paper, to point out the applications of these mathematical principles to the actual 
laws of Nature. 

A fortnight ago, I believed that no writer had ever treated of Optics on a similar plan, But 
within that period, my Tutor, the Reverend Mr. Boyton, to whom I had communicated some of 
my results, has shewn me in the College Library a beautiful Memoir of Malus on the subject, 
entitled “Traité d'Optique,” and presented to the Institute in 1807.+ 

Those who may take the trouble to compare his Memoir with mine, will perceive a difference 
in method and extent. With respect to those results which are common to both, it is proper 


to state that I had arrived at them in my own researches, before I was aware of the existence 
of his, 


Dublin, 7 S. Cumberland Street, 


December 6, 1824. 
General Principles, 
I, 
Let us conceive that from every point of a given surface 
F(a’, y, 2)=0 
proceeds a ray of light, represented by its equations 
iar dd 
mpm la PEE aaki (R) 
in which p, v are known functions of @’, y’, z' and their partial differentials = A &c., which 
we will denote by p, q, &c. 


* [See Appendix, Note 1, p. 462. This paper deals with the properties of a general rectilinear congruence, and 
with its singularities; for a recent treatment of the singularities, see Œ, Julia, Éléments de Géométrie infinité- 
simale (Paris, 1927), pp. 204-217.] | 

. t [Mémoires présentés à 0 Institut par divers savans, 2 (1811), pp. 214-302. ] 
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Such a system of rays has several properties, independent of the form of the functions F, p, v 
these properties we propose to investigate. When the equation of the surface is given, we can 
deduce from it the partial differentials p, q, &c., and therefore y, v, as functions of a’, y’, z’. This is 
a consequence of our fundamental hypothesis that from each point of the surface proceeds but 
one ray of the system. 

If it be required to find the ray which passes through a given external point æ”, y”, 2”, we 
have these three equations 

” d 


n t 
AA a d + Asc UD T 
-r r = : =y, F(a z’)=0 
N p Binat, ’ ( 5 yY , ) > 


to determine a’, y’, 2’, the point on the surface from which the ray proceeds, 
But since these equations are in general of a degree higher than unity, let the exponents of 
their respective degrees be 


M, N, ®, 
and the final equation will be of the degree marked by the product of these exponents MNỌ, 
which we will denote by n. nls, is then the number of possible solutions, and n different rays 
may intersect in one point æ”, y”, 2” 

If n be an even mambet, ai the roots may be imaginary. Then no ray passes through the 
assumed external point æ”, y”, z”. But ifn be an odd number, which requires that the equation 
of the surface should not W of an even degree, then one root at least is always real, and every 
point of space has at least one ray passing through it. 

If the three equations in a’, y’, z’ reduce themselves to two, they express not a finite number 
of points but a curve on the given surface, all the rays from which intersect in the point x’, y”, 2”, 
and compose a cone, having that point for its vertex. 

Finally if the three equations reduce themselves to one, all the rays intersect in the point 
æ", Y”, 2", which is then the focus of the system. . 

If we trace any curve, plane or of double curvature, on the given surface, the rays proceeding 
from that curve will compose a surface of rays. We will obtain its equation by eliminating 
ax’, y', 2’, between the equations of the ray (R) and those of the curve 

F(a, y ey 0, Y (a) 9,2) = 0, 
This curve I call the base of the surface of rays. There are an infinite number of such surfaces, 
corresponding to the infinite number of possible bases: but they have all one common property, 
and are represented by one common equation in partial differentials of the first order, 
dz dz 
bat v 7 1=0, (P) 

in which we are to write, instead of p, v, their values as functions of a, y, z, deduced from the 
equations of the ray and of the given surface. 

(P) has been obtained by elimination of the arbitrary function Y. It expresses that the 
tangent plane to the surface of rays entirely contains the ray passing through its point of contact. 

This equation (P) furnishes a criterion whereby to ascertain whether a proposed surface 
V = 0 is composed of rays: and if not, to determine the curve of contact, in the whole extent of 
which it is enveloped by a system of rays. For if V = 0 represent a surface of rays, the equation 


dV dV dV 
pie CP ee tr =0 = 
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is identically satisfied: and if not, V=0, v=0 are the equations of the curve of contact. This 
curve separates, on the proposed opaque surface V = 0, between the part illumined, and tke part 
obscure. 

There are two particular surfaces, or rather two naps of one surface, which though not surfaces 
of rays, yet satisfy the equation of condition (P). Their equations 
i 2=0, X =0, 
without containing an arbitrary function, or even an arbitrary constant, are nevertheless integrals 
of that equation in partial differences of the first order. 

In fact they are singular primitives of that equation: and they belong to the two caustic 
surfaces, touched by all the rays. This subject will be resumed in the IIId Article of this paper. 

To find the surface of rays passing through a given curve V=0, v = 0, we will first find its 
base, V (a’, y’, z’)= 0, by eliminating a, y, z between the equations of the given curve, and those 
of the ray. But it is necessary to observe, that since through each point of the given curve may 
pass n rays, the problem has in general n solutions: or rather the surface of rays has in general 
n branches, which all intersect in that curve. When V=0, v=0 are the equations of a given 
ray, instead of a curve, the surface of rays has only n — 1 branches: for one evidently in this case 
coincides with that given ray. 

To find the surface of rays, enveloping a given surface V = 0, we are first to find the curve of 
contact, and so reduce this problem to the last. 

These two problems contain the complete solution of the problem of shadows, The caustic 
surfaces have no determined shadow. Every other curve or surface may in general have n shadows. 

Sometimes all the rays of the system intersect. the curve V=0, v= 0. In this case, that curve 
takes the place of a caustic surface: and the equation 

V (a’,y',2)=0, 

found by elimination of æ, y, z between the equations of the curve and those of the ray, is iden- 
tically satisfied by that of the given surface F (a’,y’, z’)=0. 

The subject of caustic surfaces will be more fully treated in Article III. 


II. 


In the general equations of the ray, w, v are known functions of «’, y’, 2’ and their partial 
differentials. But these latter are themselves functions of a’, y’, 2’, derived from the equation of 
the given surface 

F(@,y, 2)=0: 
so that, all reductions made, we may denote the differentials of 2’, p, v as follows, 
dz =pda' + qdy', ` du=€da' +ndy', dv=C/da' +7'dy’, 
P, 4% & n, &, n being all functions of a’, y’, 2’. 


This being laid down, suppose that from a given ray we wish to pass to a consecutive ray 
which intersects it, we must of course differentiate its equations, making 2’, y’ vary, and z’, p, v, | 
as functions of these, but supposing dæ =0, dy=0,dz=0. That is, we must combine with the 
equations of the ray, namely 


wage Sabo ® 
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their derived equations 


p.d? =d _ v.d? —dy' _ 
BEN OA Arua M a = dp, Thanet aa = dv. (1) 


These four equations determine not only the coordinates of intersection, but the direction by 
which we must pass on the given surface, to obtain a consecutive ray intersecting the given ray. 
There are two such directions: and their properties will form the subject of the present Article. 


Eliminating z — z’ between the two last equations we obtain the following, 

(udz’ —da') dv = (vdz' — dy’) dp, (2) 
which when we substitute for dz’, dp, dv their values, becomes quadratic in w , and determines 
the directions. Before making this substitution, I will observe that when the rays are supposed 
perpendicular to the surface, on which hypothesis 

A ee a 
(pdz' + da’) dq =(qdz' + dy’) dp, 
a known form for the equation of the curves of greatest and least curvature. 


By a reasoning precisely analogous to that which is usually employed with respect to these 
curves of curvature, we infer that equation (2) is the differential equation belonging to a class of 
curves, two of which pass through every point on the given surface, and the rays proceeding from 
which intersect consecutively. These curves we will denote by (s), (s’). 

Let us now develope (2), which determines their directions. It takes the form 


this last equation becomes 


Ady”? + Bdy' da' — Oda"? =0, (3) 
in which 
A= nq +7(1—vq), 
B= O' uq + (1 — vg) —7' (1— up) — nyp, 
C= 0 (1— up) + Sop. 
dy! gee ee ad 
Suppose the roots of (8) are dg € and Pde 6”: which gives us 
A at A 
6+6 wet G 66 w 
Then the tangent to the curve of double curvature (s) has for its equations 
ham E., A WAS, Em A. A 
oa at P h phs, enp Pag Et i (t) 
and in the same manner the tangent to the other curve (s’) is given by the equations 
YY 9 8a ee j A 
Say me y eee Pe ga (t’) 


Representing by @ the angle under which the curves (s), (s’), or their tangents (t), (t’), cut 


each other, we easily deduce 
(1+ p? +97) (€—6’? = N? tan? 6, 
N being 
1+ 66" +(p+q6)(p+q6’)=N. 
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If we observe that these formule contain only symmetric functions of 6, 6’, the roots of (8), 
we will find no difficulty in calculating them in terms of A, B, O, the coefficients of that equation. 
Thus we obtain 

(1+ p? +9?) (B? +440) =(AN)* tan? 0, (4) 
AN =(1+p?) A —pqB-(1+@) 0. 

By means of this equation (4), if the law and the surface be given, we can find the angle 9, 
made by the two directions at any given point on the surface: or that angle being given, we can 
find a curve on the surface, the locus of the points at which the two directions cut under that 
given angle. 

If the surface be given but not the law, (4) is the equation of condition which that law must 
satisfy, in order that at every point of the given surface, the two directions should cut under the 
required angle 9. 

Finally if the law and the angle are given, but not the surface, (4) is the equation in partial 
differences of the latter. 

Two particular cases deserve notice: Ist. 0 =90°, 2nd. 0=0. When 0 =90°, that is, when the 
two directions cut at right angles, 

tand=o, or (1+p*?)A—pgB—(1+q)C=0. 
This condition becomes, when we substitute for A, B, C their values, 


E (L +g- up) + tp (v +g)=n'g (u +p) +n (1 +p* — vg), 
and it is satisfied of itself when the rays are perpendicular to the surface: for then 
y+q=0, p+p=0, 1+q?-pp=1+p?— rq, ==- 
When 0 =0, that is, when the two directions coincide, 
tanĝ=0, or B?+4AC=0. 
The curve on the given surface which has for its equations the two following 
F(a',y',z)=0, B+440=0, (A) 


is a limit of a remarkable nature. Every ray proceeding from a point on one side of that limit, 
is intersected by two consecutive rays: every ray from a point upon it is intersected by but one- 
consecutive ray: and finally every ray from a point on the other side of it, is not intersected by 
any consecutive ray. 


For if B? + 4A0 >0, the roots of Ady’? + Bdy' da’ — Cda'*=0 are real, and tan @ is real : 


W EF GR BO os coca ieia a aa o pete aes equal, and tan 0 = 0: 
SF Pe a nach ain bettbn tina OTTE PEI EOE NR imaginary, and tan ĝ is 
imaginary. 


This limit (A) is in general touched by all the curves (s), (s’), and its equations form a singular 
- primitive of their differential equation (3).* 


* [This is not necessarily true. (A) is a locus of cusps of (s), (s’), but is not, in general, touched by these 
curves; (A) does not, in general, satisfy (3); cf. E. Picard, Traité d Analyse, 111 (1908), chap. m1; E. Goursat, 
Mathematical Analysis (Hedrick-Dunkel), 11 (2), p. 198 et seg. The error in Hamilton’s reasoning is explained in 
the Appendix, Note 22, p. 504.] 

HMP . 45 
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But we must observe that if any one of these hypotheses respecting the values of B? + 4AC be 
for any particular surface or law impossible, the limit (A) in that case does not exist at all, or at 
least not in the sense we have described. 

If, for example, B?+ 4AC be essentially positive, or essentially negative, in the one case every 
ray is intersected by two consecutive rays, in the other case no two consecutive rays intersect. 

Again if B?+ 4AC0=0 is possible, but B?+4AC<0 impossible, the two directions coincide 
at the limit, but they do not become imaginary at either side of it. In this case the limit (A) is 
not touched by the curves (s), (s’), but cuts them obliquely. I will mention two cases of this, which 
are both of them found in Nature. 

Ist Case. When B?+4AC is an exact square: in this case the differential equation resolves 
itself into two rational factors. 

2nd Case. When B?+4AC=0 can only be satisfied by putting separately * 


A=0, B=0, C=0: 
and when these suppositions are possible. 

In the VIth Article of this Essay + will be found investigated a class of laws, including normals, 
reflected, and refracted rays, which come under this Second Case. The consideration of the point 
to which, in this case, the limit (A) generally reduces itself, is important in the Theory of Optics: 
and in the Second Part of this Essay + it will be introduced with all necessary development. 

Let us conceive the system of equations 

F(a’, y', 2')=0, Ady”? + Bdy'da' — Cda'*=0 
actually integrated. From their nature, the integral may be put under the form} 
F(a’, y', 2)=0, (c-—oP=%, 
the first being as before the equation of the given surface, and the second determining upon it 
the curves (s), (s’). $, y represent rational functions of æ’, y’, 2’, and c is the arbitrary constant. 

To complete the integral, that is, to determine the constant by the condition of the curves 

(s), (s’) passing through a given point of the surface, we have 


d =[p] +y, ce’ =[$]-[vE 
[$], [Y] representing what $, yr become at the giveu point, and c’, c” being the corresponding 
values of c: so that the equations of the two curves passing through that point are 
('-gP=h (s) and (d -pf=ẹp (s). 
These curves are distinct, coincident, or imaginary, according as 
[y]>0, =0, or <9, 
conditions which correspond to those before given, namely B? + 4AC >0, =0, or <0. 


They are altogether independent in their nature, if the function w is an exact square, in 
which case B?+4AC is an exact square also. There exists then no one curve touched by them 
all, but all of one kind may intersect in one point. 


* [This happens when the rays form a normal congruence. | 
+ [Not included in the MS.] 
t [See Appendix, Note 22, p. 504.] 
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We shall get an idea of these circumstances, if we recollect that the curves of greatest and 
least curvature on a surface of revolution are the meridian and the parallel: the latter is always 
a circle, but the former may be any plane curve; there exists no one curve touched by all of 
either species, but all the meridians intersect in the vertex of the surface. 


In general, it may happen that through some particular point of the given surface pass more 

than two curves (s), (s’): in this case the quadratic equation 
Ady’? + Bdy'da’ — Cdz'*=0 

is not proper to determine their directions: analysis shews this by the coefficients vanishing. 

If there be three directions, they are given by the equation of the third degree 

dA .dy'?+dB.dy' .da’' — dC. da? = 0, 
and in general if there be a finite number of directions = n, they are given by 
d'= A dy? +d" B .dy' .da' —d"* 0 .da'*=0, 

which becomes, at this point, an equation of the first order, and nth degree. 


The surfaces of rays proceeding from the curves (s), (s’) are all developable, because composed 
of lines which intersect consecutively. Their equations are deduced from the equations of the 
curves (s), (s’), by the method of the Ist Article. 


In the next Article, we will treat of their aretes:* in the present we will confine ourselves to 
their tangent planes, 

There are in general two such developable surfaces passing through every ray: their tangent 
planes contain the elements of the two curves (s), (s’), that is to say, their tangents (t), (t’). Hence 
the equations of those planes are 

Ist.  (w— a’) {(1—g) 6 — vp} + (y—y’) {ug —(1—pp)} +(2-7){v-u6}=0, (T) 
Ind. (#—a’) {(1—vg) 6’ — vp} + (y—y') (ug? —(1—pp)} + @— 2) {y—w6}=0; (T) 
€, € are the roots of the equation Aé*+ B6—C=0: so that 6+6'=-4, e6'=-©. The 
trigonometric tangent of the angle w under which these two developable surfaces, or their tangent 
planes, cut each other, is given by the equation 
(1— pp — vg) (1 + p? + v*) (B? +440) = N” tan? w, 
N’ being AA’ + BB’ — 00”, and 
A’=(1—pp?+v?(1+p*), 
B' =(utp)(v+q)—pqg( +e? +r, 
i O =(1—vq)? +p? (1 + @?). 

This formula is susceptible of an important reduction. In fact, if we develope the quantity 

here denoted by W’ we find it is of the form 


N” being oia z wai 

p+) E+ A- up+i)n— (1 -vg+ p) E — v (n+p): 

so that, setting aside the particular case in which 1— up — vq =0, our angular formula becomes — 
(1 + p? + v*) (B?+4+4A0)=N"" tan? o. (5) 


* (Throughout the MS. this word is spelt without the circumflex accent. ] 
45-2 
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It is easy to apply to this equation (5) the remarks that were made in page 349, on the uses 
of equation (4), The same particular cases, here also, deserve notice, namely w =0 and œ = 90°. 

When o=0, tana=0, B?+4AC=0: that is to say, the tangent planes to the developable 
surfaces coincide at the limit (A). 

When w= 90°, that is, when the developable surfaces cut at right angles, we have tan w =, 
Dp 0, 

Bvt Qo+l—pptr*)n=(l—vgt p) Holp). (6) 
This condition is always satisfied when the rays are perpendicular to the surface: for then 
v+q=0, wt+p=0, l—-ppt+=1-19¢4+ p23, n=- t. 

It is also satisfied by reflected and refracted rays. 

I will conclude this Article by shewing how some of the preceding formulæ may be simplified, 
1st, by a choice of the axes: 2nd, of the independent variables. 

lst. When our object is not to examine the general properties of the system, but circum- 
stances connected with some particular ray, we may assume that ray for the axis of z, which gives 
a =0, y' =0, p =0,v=0. On this hypothesis the equations of the ray become (R)......c=0, y=0 
their derived equations (1) become 


the condition of intersection (2) becomes 
dv. da’ = dw. dy’; 
its developed equation (3) becomes 
ndy? +(&—1') dy’ da’ — o' da’? = 0, 


the coefficients being 
A=y, B= f-7, C=¢'; 


the tangent planes become 
Ist. y—6x=0 (T), IInd. y- 6's =0 (T); 
the formula for their mutual inclination becomes 
B? + 44C = (4 — CF tan? a, 
(E-n) + 4ng = (n — Ẹ'Y tan? o, 
and the condition for their rectangularity, A = O or n=€’. 


or 


, 
In this position of the axes, ni is the trigonometric tangent of the angle, which a plane con- 
1 
taining the ray and cutting the surface in a direction corresponding to 2 makes with the 


vertical plane zæ, which latter also contains the ray. We will denote this angle by w, and we 
have, for the two tangent planes 
n tan? w + (€—7’) tana — ¢' = 0, 
which may be thus written 
n— ¢' +(€—7')sin 20 = (n + £') cos 2a, 


under which form we shall have occasion to employ it in the IVth Article.* 


* [Not included in the MS.] 
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Finally if we denote — 1, the reciprocal of the focal distance taken negatively, by A, we 


have from equation (1) 
A=€+%76=€+7 tan w. 


2nd. The other simplification, which however we shall not make use of, consists in supposing 
the independent variables to be a’ — yz’, y’ — vz’, the coordinates of the intersection of the ray 
with the horizontal plane. Denoting these by a, b, the equations (R) of the ray become 


e-a=pe, y—b=vze: 
the condition of two consecutive rays being in the same plane is 
du . db = dv . da: 
and the condition of the two heen ike right angles is 


d 
py. Het) Haat ay? + py. ot 


LEL 
The general equations of the ray being 


ol ae yy _ 
Pas ies BS (R) 


we shewed that to obtain its intersection with a consecutive ray of the same kind, we must 
combine with (R) the derived equations 


nde -W = pda + ndy =dp, PEH =o dal +n! dyf =v: 
ray ate y 


Z 


equations which we will now write thus 


æ=0 +ul, yay +T, cz=2' +T, (R) 
_ pdz? —da! _ vdz' — dy’ , nh ugh 
T= yee mra yr peer (1), I being z— 2’. 
The equation of condition 
(udz’ — da’) dv =(vdz’ — dy’) dp (2) 


was shewn to indicate two directions by one or other of which we must pass on the given 
surface, in order that two consecutive rays should intersect. Being developed it was put under 
the form 

A dy’? + Bdy'da' — Cda? = (3) 
and shewn to be the differential equation of a class of curves (s), (s’), analogous to those of 
greatest and least curvature, two of which pass through every point of the given surface, (except 
those points which are at or beyond the limit (A),) and of which we represented the integral 


equation by i 
| (p-of ay. 


We calculated the tangents to these curves, and the angle under which they cut each other: we 
shewed that the rays proceeding from them intersect consecutively, and so compose developable 
surfaces (D), (D’): we calculated the tangent planes to these surfaces, and the angle œ under 
which they cut each other. 
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These were the principal subjects of the last Article. In the present we propose to consider 
the aretes of those developable surfaces; the coordinates, loci, and other circumstances connected 
with the intersections of consecutive rays. 

To find the coordinates æ, y, z of the point in which the given ray from the point a’, y’, 2’ is 
intersected by a consecutive ray from the point 


æ' +da’, y +6da', 2'+(pt+qé) dz’, 
a’ + da’, y' + 6'da’, z! +(p+q6’) dz’, 
6, 6’ being roots of the equation 46? + B6—C=0, we have the following equations : 
æ=0 +p, ysy +r, z= +T, (R) 
1- p+iT + (nl —pq)6=0 (i), 1 pap t bT + (aT — pa) E= 0 (7). 
The three equations marked (R) are common to both points of intersection : with them we 
must combine (i) for the first point, and (i’) for the second. 


The focal distance of the surface in each direction, that is, the distance of the point of inter- 
section from the surface, measured on the ray, has for its expression. 

£{@—ay+(y—y'P+ e-t NIFA. 

To find the arete (a) of the developable surface (D), that is, the locus of the intersections of 
the rays from the curve (s), eliminate a’, y’, 2’, T between the equations of that curve (s), and the 
four equations marked (R), (i): in the same manner to find the arete (a’) of the developable 
surface (D’), eliminate a’, y’, 2’, T between the six equations (s’), (R), (i’). 

These aretes (a), (a’) are the caustics of the curves on the given surface (s), (s’): for they are 
touched by all the rays from those curves (s), (s’). These latter are the only curves on the given 
surface which have caustics. 

The tangent planes (T), (T’) of which we found the equations in page 351, are the osculating 
planes of the caustics (a), (a’): for the tangent plane to a developable surface is always the 
osculating plane of its arete. 

Every ray of the system, within the limit (A), is the intersection of two developable surfaces, 
and the common tangent of two caustics, 

The locus of the caustics (a) is a curve surface (C), touched by all the rays, and containing 
upon it all their intersections of the first kind: we will call it the first caustic surface. Its 
equation is obtained by eliminating 2’, y’, 2’, T between the equations (R), (i) and that of the 
given surface, F(a’, y’, 2’) = 0. 

In the same manner, the locus of the caustics (a’) is a curve surface (C’), touched also by all 
the rays, and containing upon it all the intersections of the second kind. We will call it the 
second caustic surface, and obtain its equation by eliminating a’, y’, 2’, between (R), (i) and 
F(¢,y',2)=0. , 

These two caustic surfaces are only, in general, two different branches of one and the same 
great surface: for if we represent their equations by }=0, >’ =0, we find that these are in 
general irrational factors of the equation 2>’ = 0. 


or from the point 


e t 
This collective equation is immediately obtained by eliminating 2’, y’, 2’, 4 between the 


equations of the surface and ray, and their derived equations : or, in our present notation, by 
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eliminating a’, y , 2’, I between the equation of the given surface, and these four other equations 
vag +r, y=y trl, z= +T, (R) 
(9D +1 — vg) ($T +1 — up) = (nT — ug) (Ẹ'T — vp). (1) 
The last equation may be thus written 
A”'T? + B’T' +0” =0, 


in which 
A” =9'F—nf', 
B” =(1—pp—vq) (+n) + (pn +4’) + we (ph + 98"), 
0” =1—pp— vg. 


We cannot give a complete discussion of this equation, until we have shewn how to calculate 
the partial differentials of the caustic surface. However some obvious remarks may be inserted 
here. 

a ! ; ; i 
If qr > 0, the roots, if real, are of the same sign, and the two foci are at the same side of the 


given surface. If T < 0, the roots are necessarily real, and of opposite signs: then the two foci 


+ 


are at opposite sides of the given surface. 

If B” =0, the roots, if real, are equal in quantity, but different in sign: then the two foci are 
at opposite sides, and at equal distances. 

if B”2— 4A” 0” =0, the roots are equal, and the foci coincide. This condition being 
developed is found to be the same as B? + 440 = 0. 

If A” =0, one of the roots is infinite, and two consecutive rays parallel. In this case 
dp = 0, dv =0: that is 

Eda'+ndy' =0, ¢'da' +n'dy'=0: 

either of these equations determines the direction on the surface in which one of its focal 
distances is infinite. 

Every ray from the curve on the given surface, determined by the equation 

A” =1'F— n% =0, 


is asymptot to two branches of a caustic situated on different sides of the given surface : and the 
surface of rays from the same curve is asymptotic to two branches of a caustic surface, situated 
likewise at different sides of the given surface. If A” =0 be identically satisfied, one of the 
developable surfaces is cylindrical. 

Let us now consider the problem “To find the partial differentials of the caustic surface,” its 
equation in #, y, z being given as the result of elimination of a’, y’, z’, IT between the five equations, 
F(a’, y', 2) =0, (R) and (I). 

Instead of stopping at this, I will propose at once the more general problem “Given a system 
of m equations between m +n vatiables, to find the partial differential coefficient of any one of 
them, relative to the variation of any other.” 


Let us represent the coefficient sought by M : and conceive m —1 variables eliminated 


between the m equations, so as to leave one equation between n + 1 variables, including w and z. 
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Of these we will consider z as the dependent variable, and the n others as independent. If now 
we differentiate, suppressing all the differentials except those of # and z, the resulting equation 
in dz, dz gives the coefficient. 

But the elimination which we have here supposed, may be impossible, and is never necessary. 
In fact if we attend to the principle of the method, we find that it consists in supposing constant 
n— 1 out of the n independent variables: on which hypothesis the ratio dz: dw is the coefficient 
sought. 

If then we differentiate the m given equations, supposing n — 1 of the variables constant, the 
m resulting equations, being linear equations between m partial differentials, contain an easy 
solution of the problem. 

Perhaps I may have not succeeded in expressing myself with clearness, on a subject so very 
general; but I hope that the preceding remarks will be understood by observing their application 
to the case before us, 

In this case, we have a system of five equations between seven variables, namely a, y, z,#’, y’,2’,T': 
consequently two are independent; for these we will take a, y. Now differentiating, and supposing 


dz 
y constant, we obtain five linear equations between As dp da dea? five partial differentials 


of which we only want the first. But since these are all linear equations there is no difficulty, 
except the length of the calculation, (which in this case is abridged by some reductions that 


present themselves,) in finding each separately. In a similar manner I have calculated A and 
these are the results, after all reductions :* 


Rikat Trpia Tika g 
aQ Wi- M-a 


PEL iS Mt Sell 2 eel 
vP CT-—yp Tl- 


P, Q here denote the partial differentials s ; T , of the caustic surface. 


The two expressions thus found for an l are equal among themselves, in consequence of 


equation (I): and it is easy to see that they are the reciprocals of those for i . P 4 d 


* [These results may easily be obtained from the equations 
a=a+pl, yay +T, 2=7 +T. 
Taking an arbitrary displacement dz’, dy’, dz’ on the surface F(x’, y’, 2’) =0, we find 
dx—pde_ (1-pP)dx-pQdy _ (1+{T—mp) da’ + (nT — pg) dy 


dy—vdz —vPdx+(1—vQ) dy — ~ (T= vp) da’ +(1+n'T— vq) dy 
But, since T satisfies (I) (p. 355), this fraction is independent of dy’/dz’. Therefore 


1-pP_ —n@ _14+(T—-pp_ _ pg ] 
—vP 1-1Q (T—vp ~ T4y'T-vq" 
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Hence we deduce 


1Q—1 pP- 
Po TE a 
that is, 
1— uP —rQ=0. 


Now this equation in partial differences of the first order, belonging to the caustic surfaces, is 
precisely the same as that which we found in the first Article, to belong to all surfaces of rays, 
by elimination of an arbitrary function. 

And it is on the analysis here described, that I founded the assertion in that Article, that the 
equations of the two caustic surfaces, which we will still consider as distinct, and represent by 
X =0, >’ =0, are singular primitives of the equation in partial differences 


dz dz 
wer haha ag 


We have already considered these surfaces as formed by two different generations, and given 
two corresponding methods to find their equations. 

In the first method we obtained each separately, considering it as the locus of the caustic 
curves (a) or (a’) of its own kind. In the second method we obtained their equations collectively, 
considering the great surface of which they are in general different branches, as the one locus of 
all the intersections of consecutive rays. 

In each method, their equations, whether separately or combined, are presented as the result 
of elimination. 

Now it is a received principle in Analysis, that whatever is represented by a system of equa- 
tions, may be represented by an infinite number of systems equivalent: in the same manner as, 
in Mechanics, forces which have no resultant, may be replaced by an infinite number of equivalent 
systems. — 

Conformably to this general principle, I am going to shew that the caustic surfaces may have 
an infinite number of different generations, among which we are at liberty to choose, in any 
particular case, that which may offer the greatest facility. 

Conceive any curve whatever, with an arbitrary parameter, traced on the given luminous 
surface F(a’, y’, z’)=0: the rays proceeding from that curve will compose a surface of rays with 
an arbitrary parameter, which we will represent by 


D (x, Y, z, €) =0 =u, 
If we consider, a consecutive surface of the same kind, 


du 
ut+—7 .de=0, 


we find that it intersects the former in three different ways. 
Ist. In the ray proceeding from the intersection of their consecutive bases. 


2nd. Every curve of the form (s) cuts the two bases in two points: the consecutive rays, one 
on each surface, proceeding from these two points, intersect each other: the locus of such inter- 
sections is a curve common to the two consecutive surfaces of rays, and this curve is on the first 
caustic surface. 


HMP 3 46 
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8rd. By a similar reasoning, the two consecutive surfaces of rays have another curve of 
intersection on the second caustic surface. 
The ray and the two curves have for their common expression 


and if we eliminate e between these, we obtain in a, y, z an equation 
E2, S =O: 


3 =0, X} =0 are the caustic surfaces; H =0 is a surface of rays enveloping all those of the 
assumed form Ẹ (x, y, z, €)= 0. 

In general every surface of rays envelopes both caustic surfaces. But there exists one remarkable 
exception which it is necessary to notice. The developable surfaces of the first kind (D) do not 
touch the first caustic surface (C): and in general the developable surfaces do not touch the 
loci of their own aretes: a circumstance sufficiently evident, because their tangent planes 
osculate to those aretes: but I believe it may be useful to shew how the preceding reasoning fails 
in this case. . 

In fact it is evident that a curve of the form (s) does not cut the bases of two consecutive 
developable surfaces (D), though it may coincide with one of them: so that if we represent by 
D=0, D'=0, the general equations of the two developable surfaces, with the arbitrary constant, 


dD 


and eliminate c between D =0, P ia 0, we obtain an equation Z>’ = 0, which does not contain 
U 

the first caustic surface. In the same manner, eliminating c between D’=0, ae =0, we 

obtain ZX = 0. 

The common factor Z=0 belongs to a remarkable surface, the principal properties of which 
I will rapidly enumerate. Ist. It envelopes all the developable surfaces of both kinds (D), (D’); 
2nd. it is a limiting surface, analogous to the limiting curve (A), since through every ray at one 
side of it pass two developable surfaces, through every ray upon it passes but one, and through 
every ray beyond it passes none; 8rd. it has for base the limiting curve (A); 4th. it envelopes 
both the caustic surfaces in one and the same curve of contact; 5th. it is itself developable and 
this curve of contact is its own arete; 6th. that arete is a curve of rebroussement on the caustic 
surfaces, in which the two branches touch each other and at which they terminate ; 7th. the same 
arete is itself a caustic curve, and the one in which the density of light is greatest; 8th. it is 
touched by all the caustic curves of both kinds, and it is the caustic of the limit (A).* 

This is a proper place to point out some developments, and some demonstrations, which may 
assist in conceiving the connexion between the curves (s), (s’) and the limit (A), and at the same 
time between the surfaces of rays proceeding from them, that is, the developable surfaces (D), (D’) 
and the limiting surface (L). 

We have represented the general form of the equations of the curves (s), (s’) by 


F(a’, y',2)=0, (¢—cP=y¥, 


* [Some of these properties of (L) (the surface of rays through (A)) are derived from. the assumption (see 
footnote to p. 349) that (A) is touched by the curves (s), (s’). In general (L) is not touched by the developables of 
the system, nor is it itself developable, but it does touch both the caustic surfaces in a single curve. ] 
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of which we need only attend to the latter, because the former belongs to the given surface. If 
we eliminate the arbitrary constant by differentiation we find the following, 


4y . dẹ’ = dẹ’, (a) 
which must be identical with the original differential equation of the curves, 
Ady + Bdy'da’ — Oda? =0; 
expanding and comparing, we find 


nam (iy) =- (37) > 


0 oe 


Leones (i) - ayh) 


dy dọ dp dp 
A (B® + 4A0) = 164. Fed ‘dy dy’ St, 
A being an indeterminate coefficient, which proves that when y > 0, =0, <0, or an exact square, 
we have at the same time B?+4AC >0, =0, < 0, or an exact square. 
It is evident that the differential equation (d) has a singular solution, y =0, dy=0. These 
hypotheses belong to the limit (A), and shew that this latter curve is touched by all of the form 
(s), (s’).* It has however only a contact of the first order with each: for if we differentiate (d), 


and put y=0, dy = 0, we find 
dp . (dp — 2d?) =0: 


the first factor belongs to the limit, the second to the curve: so that nie have at their common 
point a common tangent, but not a common osculating circle. 

Since at the limit dy = 0, d*y = 2d¢*, or d'y >0, y is a minimum, Only one curve of the 
form (s), (s’) passes through each point of the limit. In fact these two curves are in general of 
the same nature, and no otherwise differ than by the value of the constant c: which at the limit 
has but one value. 

Perhaps a simple geometrical illustration may be useful, to make these circumstances easy of 
conception. Imagine a series of circles of given radius, all touching externally a given circle, in 
its own plane: through every external point there pass in general two such circles; these we 
may conceive to represent the curves (s), (s’). Through every point upon the given circle passes 
but one, namely the circle of the given radius, which touches it externally at that point. In this 
case the given circle is the limit (A), and no curve of the kind we are considering passes through 
any point within it. It is necessary however to observe that in this case there is another limit, 
also a circle, touched by all of the given form, through every point of which passes but one, and 
through every point beyond which passes none: this is a peculiarity of the instance we have 
selected, and is not found in the general case. ` 

But when the function is an exact square which we will represent by 0%, the differential 


equation 
dnp. dẹ? — dyy* =0 (d) 
* [This is not true in general. See Appendix, Note 22, p. 504.] 
46-2 
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resolves itself into two factors: one of them is Q? = 0, the other d¢?—dQ?=0: the first belongs 
to the limit, which is still a particular solution, and the second belongs to the curves which on this 


hypothesis become 
¢+0X=c (s), @-QNQ=c (s): 


in this, as in all cases, at the limit the two values of the constant become equal, but the curves 
themselves do not touch either the limit or each other.* In this case then, (A) is no longer a 
limit, properly speaking, nor is (L) a limiting surface. 

In this case also the curves (s), (s’) and the surfaces (D), (D’) are altogether independent. 
Nor does there then exist any one surface touched by all the developable surfaces of both kinds: 
but as we observed in page 350, that in this case the curves of one kind may all pass through one 
common point, so we will now observe, that all the developable surfaces of one kind may touch 
one common surface, or pass through one common line: which surface or line is then to be 
considered as one of the caustic surfaces, in this case independent. 

It is proper to take our illustrations from the properties of normals to a surface, because these 
being treated of by all elementary writers, are familiar to the conception, and because Monge has 
taken them for the principal subject of his profound and beautiful “ Analyse.” Suppose then the 
rays all normal to a surface of revolution: one of the developable surfaces is the meridian plane, 
the other is a right cone, having its vertex in the axis: their natures are different; all those of 
the first kind pass through a common line, namely the axis: all those of the other kind touch one 
common surface, namely the one formed by the evolute of the meridian revolving round the axis, 


The angle under which the developable surfaces cut the loci of their own aretes, is the same 
as that under which they cut each other, namely œw, calculated in the IInd Article. This is also 
the angle under which the osculating planes of the caustic curves (a), (a’) cut their own caustic 
surfaces (C), (C’). Hence it follows that when w = 90°, the caustic curves are the shortest between 
two points on their caustic surfaces. This is the case for normals, for reflected and for refracted 
rays. In the Vth Article} we will prove that, in the first of these three cases, the caustic curves 
are the shortest on another kind of locus, which consequently touches the caustic surface in the 
whole extent of the caustic curve. 


Whenever the condition œ = 90° is satisfied, a thread drawn in the direction of the ray, and, 
at the point where it meets either of the caustic surfaces, freely wound round that surface, would 
trace upon it the caustic curve corresponding: and conversely such a thread freely unwound 
would trace out on the given surface the curve (s) or (s’): so that the caustics are analogous 
to evolutes, 


In fact evolutes are a particular case of caustics, when the latter, as in this paper, are 
considered in all their generality. In the case of evolutes, the rays are perpendicular to the 
surface: the curves (s), (s’) are the curves of greatest and least-curvature on the latter: the 
developable surfaces (D), (D’) are the normal surfaces, with which Monge occupied himself 
so much, and which cut at right angles: the caustic surfaces (C), (C’) are the loci of the centres 
of curvature, which answer to the foci, and finally, the caustic curves (a), (a’) are the evolutes of 


* [There is some confusion in the argument, since (unless 4 = B=C=0 on (A)) the equation 
A dy? + Bdy dx — Cdz'?=0 
(which defines (s), (s’)) is satisfied by a unique direction when B?+4AC=0. The curves (s), (s’) therefore touch 
each other at the points where they meet (A).] 
+ [Not included in the MS.] 
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the curves of greatest and least curvature, which Monge proved in a similar manner to be the 
shortest lines between two points on their own loci. 

Monge also proved that if from any point of space we could see the two surfaces which are 
the loci of the centres of curvature, their profiles would appear to cut at right angles. I will only 
stop to observe that if œ be any constant angle, the profiles of the two caustic surfaces would 
appear to cut under that angle. 

We may now resume the discussion of the equation of focal distance in page 355, 

A” T?+ B’T +0" =0, 
which was there left incomplete. We there discussed the cases A” =0 and B’=0. The 
hypothesis 
0” =0=1-—pp— rg 
belongs to a curve on the given surface, in the whole extent of which, one focal distance being 
= 0, one caustic surface intersects the given surface. 

This is a remarkable curve: let us consider some of its properties. 

Ist. It is, as we have seen, the intersection of the given surface with a caustic surface. 

2nd, But it is at the same time a curve of contact between the given surface and that 
caustic surface: for if we make 

l—pp-—vg=0, T=0 
in the expression calculated in page 356 for the partial differentials of the caustic surface, P, Q, 
we find P =p, Q=4. 

8rd. It is a singular curve on the caustic surface, for if we calculate the partial differentials, 
second order, of that surface, we find them, at that curve, in general infinite.* 


4th. The equation which in general determines tan w, disappears when 


1—pp—vq=0. 
5th. On the same hypothesis the equation which determines the directions of the curves 
(s), (s’), namely 
A dy’? + Bdy'da' — Odg? = 
resolves itself into two rational factors, of which one being ee — vdæ' = 0, that is, v — u6 = 0, 
changes the equations in page 348, for the tangent (t) to the curve (s), into these, 
y-y _v 2-2 ppt+rq_1 


B af p’ g1 g m pe’ (t) 


because in the present case 1—pp—vg= 0: and these equations (t) shew that the tangent 
to one of the curves (s), (s’) is in this case the ray itself. 


6th. Finally the curve 
F(a',y',27)=0, 1-pp—vrg=0, 
is that in which the given surface is itself enveloped by a surface of rays. 
It remains to consider the case B”? — 4A”0” = 0, 


* [Hamilton must have been led to this conclusion by an error of calculation, The partial derivatives of the 
second order for the caustic surface are infinite only if 
2p (p+ lq) + 2v (np +'q) + pr + Quvs+v2t=0.] 
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In general when we seek the intersection of the two caustic surfaces, or of the two branches 
of the great caustic surface, which is composed of them both, we suppose the two foci to coincide, 
and therefore the equation in T to take equal roots.* 

This hypothesis furnishes the condition B”? —-4A’’C’’ =0: which being developed is found to 
be the same as B?+4AC=0, Hence we conclude 

Ist. the curve of intersection is also a curve of contact: for B?4+4AC=0 gives tan w =0; 

2nd. this curve is therefore one of rebroussement on the great caustic surface ; 

3rd, the surface of rays passing through it has for base the limit (A), [and] 

4th. is therefore the developable limiting surface (L);+ 

5th. the curve itself is the arete of that surface (L) and has all the properties mentioned in 
page 358: that is to say 

6th. it is the caustic of the limit (A), and touched by all the caustics of both kinds ; 

7th. for this caustic, the density of light is a maximum: we will call it the principal caustic. 

It does not suit the plan of this paper to enter into a detail of exceptions: we will however 
mention the most extensive and important. 

Ist. When B?+4AC is an exact square. 

IInd. When B? +4AC=0 can only be satisfied by putting separately A =0, B=0, C=0, 
and when these suppositions are possible. 

These are the same exceptions that we mentioned in page 350, respecting the nature of the 
limit (A). We are going to consider them separately. 

In the Ist Case, the caustic surfaces are not branches of one and the same surface, but 
are absolutely distinct and independent. They do not touch at their intersection, nor do they 
terminate there. That intersection is not a caustic; nor is the surface of rays (L) passing through 
it, developable; nor a limiting surface: nor the base (A) a limiting curve. 

However the intersection of the two caustic surfaces is still an important curve, since for 
every point of it two foci coincide, and the density of light is greatest. 

In the IInd Case, the limit (A) generally reduces itself to a point: and there exists a class 
of laws which render not only the equation itself, 


A dy’? + Bdy'da' — Ida? = 0, 
but all its derived equations identical at that point: so that the ray from that point is intersected 


by an infinite number of consecutive rays. 
Moreover they all intersect it in one and the same point, for the equation 


A” T?+ B"T +0” =0 


takes equal roots: to this point the principal caustic in this case reduces itself: we will call it 
the principal focus. The point on the given surface to which in this case the limit (A) reduces 
itself, and from the portion of the surface near to which the light converges accurately to the 
principal focus, I will call the vertex of the surface. 


* [Hamilton leaves out of consideration an intersection for which a common point is a focus on two different 
rays.] 
+ [This, and the properties which follow, are based on the assumption that (s), (s’) touch (A).] 
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These two points exist in Nature, and it is obvious how important the consideration of them 
must be. It will of course be introduced in the second part of this paper,* whenever I have 
leisure to complete it. 


To find the coordinates of the vertex of the given surface, we may employ any three of the 
four equations 


F(a',y', 2')=0, n'ug+n(l-vg)=0, Sup +e" (l—mp)=0, mq =nvp: 
and to find the coordinates of the principal focus, a, y, z, we have the following, 


rapa, w= +T, y=y +T, z=2'+T. 


g 
I is the projection of the focal distance on the axis of z; that distance is 
———————_ 1 
2 = — 
TVI + e+= A 
I had prepared an Article on the intersections of rays at finite distances from each other, 


which was to have been inserted in this place : but the comparative unimportance of the subject, 
and the press of other matter, have induced me to relinquish that intention. 


* [Not included in the MS.] 
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